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1 Introduction 



The problem under consideration in this paper is the fohowing. 

Problem 1.1 Let T-l he a real Hilbert space, fix e S ]0, 1[, and let (/3n)neN be a sequence 
in ]0, 1 — e]. For every n G N, let : ?^ — t- be a firmly nonexpansive operator, let 
Rn'. domi2„ C Ti ^ Ti he a pseudo contraction such that (Id —Rn) is /3n~Lipschitz continu- 
ous, and let S be a closed convex subset of %. The problem is to 

find X G 5 such that (Vn € N) TnRnX = x. (1-1) 

The set of solutions to (1.1) is denoted by Z. 

As will be seen subsequently, this formulation models a broad range of problems in non- 
linear analysis. Methods can be found in the literature to solve Problem 1.1 in special cases. 
Thus, when S = H.^ Rn = Id, and Z 7^ 0, algorithms can be found in [1, 2], and when 
5 = "H, T„ = Id, and i?„ = ii, where i? is a Lipschitz pseudo contraction from a convex set 
C into itself, methods can be found in [3, 4, 5, 6]. Since the composition between a firmly 
nonexpansive operator and a Lipschitz pseudo contraction is not a pseudo contraction in 
general. Problem 1.1 can not be solved by the methods mentioned above. The purpose of the 
present paper is to provide an algorithm for solving Problem 1.1. It involves four elementary 
steps at each iteration n: the first three steps are successive computations of operators 
Tn, and Rn-, and the last step is an outer approximation of the constraint. The latter is given 
by the projection onto a half-space containing S. In Section 2 we propose our algorithm 
and we prove its weak convergence to a solution to Problem 1.1. In Section 3 we study an 
application to monotone inclusions under convex constraints, and obtain an extension of a 
result of [7]. Finally, in Section 4, we study an application to equilibrium problems with 
convex constraints. 

Notation 1.2 Throughout this paper % denotes a real Hilbert space, (• | •) denotes its scalar 
product, and || • || denotes the associated norm. For a single- valued operator R: dom.R C 
H ^ H, the set of fixed points is Fix i? = |x G 7^ | x = R is x~Lipschitz continuous 

for some x ^ ]0, +oo[, if it satisfies, for every x and y in \\Rx — Ry\\ < x\\x — y\\, R is 
pseudo contractive if it satisfies 

(Vx G domi?)(Vy G domi?) \\Rx - Ryf < \\x - + ||(Id -R)x - (Id , (1.2) 

R is firmly nonexpansive if it satisfies 

(Vx G domi?)(V2/ G domi?) \\Rx - Ryf < \\x - yf - ||(Id -R)x - (Id -R)yf, (1.3) 
or equivalently, 

(Vx G domi?)(V?/ G domi?) {x - y \ Rx - Ry) > \\Rx - Ryf , (1.4) 
and R is x-cocoercive if X-R is firmly nonexpansive. 
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2 Algorithm and convergence 



At each iteration n G N, our method for solving Problem 1.1 involves an outer approximation 
to S and separate computations of the operators T„ and Rn- Each approximation is computed 
by the projection onto a closed affine half-space containing S, and errors on the computation 
of the operators are modeled by the sequences {an)neN, {bn)neN, and (cn)neN- 

Algorithm 2.1 Let (T„)„gN, {Rn)neN, and S be as in Problem 1.1. For every n E N, let 
Qn ■ H ^ H he the projector operator onto a closed affine half-space containing S, let {an)neN, 
{bn)nen, and (c„)„gN be sequences in 7i such that ||a„|| < +oo, ^neNpuW < +oo, 



and ^„ 



< -|-oo. Moreover, let e G ]0,1[, let (A„)„gN be a sequence in [e, 1], let 



xq G domi?0) and consider the following routine. 



Un -^n-^n ~l~ 
Qn — TnUn "H 

If g„ ^ dom. Rn stop. 
Else 

— RnQn ~\~ C-n 
^n — Xn Un ~l~ 
2;^+! — Xn ~\~ ■^ni^Qn^n Xn^ 

Ifx„+i ^ domi?„+i stop. 
Else n = n + 1. 



(2.1) 



Our main result is the following. 

Theorem 2.2 Suppose that Z ^ in Problem 1.1 and that Algorithm 2.1 generates infinite 
orbits {xn)n&i o,nd {zn)nen such that 



(Vx e U) < 



Xkn ~^ X 

Xn '-^n-^nXn ^ 

^n Xn ^ 

^n Qn'^n ^ 



Then {xn)ne'N converges weakly to a solution to Problem 1.1. 



x£ Z. 



(2.2) 



Proof. Set 

(Vn G N) yn = RnXn, Qn = TnVn, and f„ = RnQn, 

fix z G Z, and let n G N. Note that, since z £ S, we have 

Z = PsZ = QnZ = TnRnZ = RnZ + (H -Rn)TnRnZ. 



(2.3) 



(2.4) 



In addition, it follows from [8, Theorem 1] that (Id —Rn) is monotone, which yields 
((Id —Rn)qn — (Id —Rn)z \qn — z) > 0. Therefore, we deduce from (2.4), (2.3), and the 
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firm nonexpansivity of r„ that 

2{qn - Z I (Id -Rn)Xn " (H -Rn)qn) = -2(5n " ^ | (W -Rn)qn " (M -Rn)z)) 

+ 2(5n - Z \ Xn - z) - 2{qn - Z \ - 

< 2(g„ - z\Xn- z) - 2{Tnyn " | TJn - RnZ) 

< 2{qn - Z \ Xn - z) - 2\\Tnyn " T„i?„z||^ 



(2(g„ - Z \ Xn - Z) - WQu - z\\'^) - WQn - Z 



2 ll~ ^ l|2 ||~ ,||2 



<||x„-2;|| -||g„-2;„|| - ||g„ - z|| . (2.5) 

Hence, since sup^gj^ < (1— e)^ < 1 — e, it follows from (2.3) and the /3n-Lipschitz continuity 
of (Id -Rn) that 

\\Xn -Vn + Tn- z\[^ = ||5n - 2 + (x„ - y„) - (g„ - f„)||^ 

= ||gn - 2 + (Id -Rn)Xn - (Id -i?„)gn||^ 

= - zf + II (Id -i?„)x„ - (Id -Rn)qnf 

+ 2(g„ - 2; I (Id -Rn)Xn - (Id -Rn)qn) 

< Mn - zf + f3l\\qn - Xnf 

+ 2(g„ - 2; I (Id -Rn)Xn - (Id -Rn)qn) 

< \\Xn - Z\\'^ - (1 - /3^)||5n - XnW"^ 

< \\xn - zf - e\\qn - Xnf , (2.6) 

which yields 

\\Xn -Vn + rn - z\\ < \\x„ - z\\. (2.7) 

We also derive from (2.1) and (2.3) the following inequalities. First, ||y„ — y„|| = ||a„||, and 
since T„ is nonexpansive, we obtain 

Ikn - qn\\ = WTnVn + K - r„y„|| < ||y„ - y„|| + ||6„|| = ||a„|| + ||6„||. (2.8) 

In turn, it follows from the /3n,-Lipschitz continuity of (Id —Rn) that 

ll^^n '"nil — ll-^n'Zn ~l~ Cjj -RnQ'nII 

< ||(Id -Rn)qn - (Id -Rn)qn\\ + |kn " qn\\ + ||c„|| 

< (l+/3n)|kn 9rt 1 1 ~l" 1 1 Cji 1 1 

<2(||a„|| + ||6„||) + ||c„||. (2.9) 

Altogether, if we set 

en = yn-yn + r„ - f„, (2.10) 

we have 

llCnll = ||yn -Vn + Tn- TnW < WVn " ^nH + l^n " J^nH < 3||a„|| + 2||6„|| + ||c„||, (2.11) 



4 



and therefore X^fceNll^fcll < +00. Hence, from (2.1), (2.4), the nonexpansivity of Qn, and 
(2.7) we get 

||Xn+i - z\\ = 11(1 - \n){Xn - z) + \n{QnZn " QnZ)\\ 

< (1 - A„)||x„ - z\\ + XnWQnZn " QnZ\\ 

< (1 - \n)\\Xn - z\\ + \n\\Zn " ^^H 

< (1 - A„)||x„ - z\\ + A„(||3;„ -Vn + rn- z\\ + ||en||) 

< ||x„-z|| + ||e„||, (2.12) 

and we conclude from [9, Lemma 3.1] that 

= sup \\xk — z\\ < +00. (2.13) 

fceN 

Thus, from the convexity of || • |p, the firm nonexpansivity of Qn, (2.4), (2.1), (2.10), (2.6), 
and (2.7) we have 

- zf < (1 - Xn)\\Xn - zf + XnWQnZn " Qnzf 

< (1 - A„)||x„ - z||^ + Xn{\\Zn - - \\Zn " Qn^nH^) 

< (1 - A„)||x„ - zW"^ + A„(||x„ - Vn + rn- zW"^ + ||e„||^ 

+ 2 \\Xn Vn ~l~ ^n z\\ ||Cn|| || Zn Qn Zn \ \ ) 

< (1 - Xn)\\xn - zW"^ + A„(||x„ - z||^ - e\\q„ - a;„|p 

+ l|en|r + 2|| Xn z\\ ||Cn|| \\Zn QnZnW ) 
— \\Xn ■^ll £ WQu Xn\\ £ ll^n Qn^nll + Vm (2-14) 

where rjn = ||en|P + 2C||era|| satisfies X^fceN'?^ < +00. Hence, from [9, Lemma 3.1] we deduce 
that 

'^WTkRkXk- XkW'^ = '^\\qk- XkW'^ < +00 and ^ H^fc - Qfc^fclP < +00, (2.15) 
fceN fceN feeN 

and therefore TnRnXn — Xn = Qn — Xn ^ and Zn — QnZn 0. Thus, it follows from (2.1) 
and the nonexpansivity of r„ that 

1 1 Zn Xn II — II ^n Vn \ \ 

— 1 1 ^n Un ~l~ (^n \ \ 

— Il^n^n 2^n2Jn|| ~l~ ll^nH 

— 1 1 Qn Xn 1 1 ~1~ 1 1 671 1 1 

^ 0. (2.16) 

Altogether, since (2.2) asserts that all the weak limits of the sequence {xk)keN are in Z, the 
result follows from [9, Theorem 3.8]. □ 



3 Monotone inclusions with convex constraints 

We consider the problem 

find X e S such that e Ax + Bx, (3.1) 
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where A: % 2^ and B: domS C % ^ % are maximally monotone, and S C H is 
nonempty, closed, and convex. When B is cocoercive, dom5 = T-L, and S = T-L, (3.1) models 
wide variety of problems in nonlinear analysis, and it can be solved by the forward-backward 
splitting method [10, 11, 12, 13, 14]. However, in several applications these assumptions are 
very restrictive. If the cocoercivity of B is relaxed to Lipschitz continuity, (3.1) can be solved 
by the modified forward-backward splitting in [7]. We propose an extension of this method 
for solving (3.1) with a finite number of convex constraints. In addition, our method allows 
for errors in the computations of the operators involved. 

Notation 3.1 For a set-valued operator A: T-L — )• 2^, domA = {x G | Ax ^ 0^ 
is the domain of A, zer A = |x G ?^ | G ^x} is its set of zeros, and graA = 
[{x,u) ^ T-L X Ti I n G Ax^ is its graph. The operator A is monotone if it satisfies, for 
every (x, u) and (y, v) in gra j4, {x — y \ u — v) > 0, and it is maximally monotone if its graph 
is not properly contained in the graph of any other monotone operator acting on H. In 
this case, the resolvent of A, J a = (Id +A)~'^, is well defined, single- valued, dom = T-L, 
and it is firmly nonexpansive. For every a G M, the lower level set at height a of a func- 
tion f : T-L ^ ]— oo,-|-oo] is the closed convex set lev<c / = {x G | f{x) < q} and the 
subdifferential of / is the operator 

df:'H^2'^:x^{u€'H\{yy€'H) {y - x \ u) + f{x) < f{y)} . (3.2) 

Now let C be a nonempty subset of Ti. Then int C is the interior of C and if C is nonempty, 
convex, and closed, then Pq denotes the projector operator onto C, which, for every x £ 7i 
satisfies ||x — Pcx\\ = min^gc" \\x — y\\ = dc{x), where dc denotes the distance function of C. 
For further background in monotone operator theory and convex analysis see [15]. 

Problem 3.2 Let A: T-L — )• 2^ and B: domi? C Ti ^ Ti he two maximally monotone 
operators such that dom A C dom B and suppose that A + B is maximally monotone (see 
[15, Corollary 24.4] for some sufficient conditions). For every i G {1, . . . ,m}, let ?^ — t- M 
be lower semicontinuous and convex, denote by = lev<o /iH - • •nlev<o fm 7^ 0, and assume 
that S C domB and that B is ^-Lipschitz continuous on 5 U dom A, for some x ^ ]0, +oo[. 
The problem is to 

X G zer(^ + B) 

(3.3) 

l/m(x) <0. 



find X G7i such that < 



Problem 3.2 models various applications to economics, traffic theory, Nash equilibrium 
problems, and network equilibrium problems among others (see [16, 17, 18] and the references 
therein). 

In the particular case when m = 1, fi = dc, and C C 7^ is a nonempty closed convex 
set, an algorithm for solving Problem 3.2 is proposed in [7], without considering errors in 
the computations and assuming that Pc is easily computable (see also [19] for an approach 
using enlargements of maximally monotone operators). However, since Ps is not computable 
in general. Problem 3.2 can not be solved by this method. We propose an algorithm for 
solving Problem 3.2 in which the constraints /i < 0,. . ., /m < are activated independently 
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and linearized, and where errors in the computation of the operators involved are permitted. 
For the implementation of this method we use the subgradient projector with respect to 
/ G ro(^), which is defined by 



X 



X, 



fix) 



\u\\ 



if/(x) >0; 
otherwise. 



(3.4) 



where u G df{x), and the function i : N — )• {1, . . . , m} : n i— )• l+rem(n— 1, m), where rem(-, m) 
is the remainder function of division by m. 

Algorithm 3.3 For every i G {1, . . . , m}, denote hy Gi: % ^ % the subgradient projector 
with respect to /j. Let {ei^n)n&]-, {^2,n)n&i-, and {e^^n)n&'^ be sequences in % such that 

II < +00, and ||e3,„|| < +oo. Let e G ]0, l/(x + l)[, 

let (7n)nGN be a sequence in [s, (1 — let xq G domB, and let (xn)nGN be the sequence 

generated by the following routine. 



(Vn G N) 



yn = Xn- ln{BXn + ei,„) 

Qn = J-y^AiUn + e2,n) 

rn = Qn- ln{Bqn + e3,n) 



(3.5) 



Remark 3.4 In Algorithm 3.3, the sequences {ei^n)n&i and {e^,n)n&i represent errors in the 
computation of the operator B. In addition, we suppose that the resolvents {J'y„A)n<m can 
be computed approximatively by solving, for every n G N, the perturbed inclusion 



find q € H such that Vn — Q + ^2,n G jnAq. 
Proposition 3.5 Suppose that 



(3.6) 



IJranGjCdomB and S nzeT{A + B) ^ 0. 



(3.7) 



1=1 



Then Algorithm 3.3 generates an infinite orbit (x„)„gN which converges weakly to a solution 
to Problem 3.2. 



Proof. Set 

(VnGN) /3„ = 7„x, Tn = J-,„A, and i?„ = Id -7„5. (3.8) 

Note that (/3„)„gN is a sequence in ]0, 1 — e] and, for every n G N, T„ is firmly nonexpansive 
and Id —Rn = "YnB is /3„-Lipschitz-continuous and monotone. Hence, it follows from [8, 
Theorem 1] that the operators {Rn)neN are pseudo contractive. In addition, note that x G 
zer{A + B) 4^ (Wn G N) x - jnBx G x + 'jnAx <^4> (Vn G N) x G Fix T„ii„. Altogether, 
we deduce that Problem 3.2 is a particular case of Problem 1.1 and 

Z = Snf~]FixTnRn = Snzer{A + B)^0. (3.9) 

neN 
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Now let us prove that Algorithm 3.3 is a particular case of Algorithm 2.1. Set 

«n = -7nei,n r ^ _ ^ 

(Vn G N) {bn = J'r„Aiyn + e2,„) - J^„Ayn and < J" _ (3-10) 

_ \^Qn — (-'i(n)- 

Then, since sup„gi^7„ < we have XlngNll'^"!! < and ZlneN ll'^nll < +°*^- More- 

over, from the nonexpansivity of (<-^'y^A)nGNi 

we deduce that X^neNll^^ll < +00, and, 
for every x G Ti and n G N, Qn^; is the projection onto the closed affine half-space 
{y G T-L I {x — y \ u) > for some u G <9/i(„)(x), which contains lev<o /i(ri) D S". 

On the other hand, xq G domi? and since, for every i G {1, . . . ,m}, ranGj C domi?, it fol- 
lows from (3.5) that, for every n G N \ {0}, Xn G domS. In addition, g„ = Jj„A{yn + e2,n) £ 
domvl C domi?. Altogether, from (3.8) and (3.10), we deduce that Algorithm 3.3 is a 
particular case of Algorithm 2.1 and that it generates an infinite orbit {xn)nm- 

Let us prove that condition (2.2) holds. Suppose that Xk„ x, x^ — T^RnXn — ^ 0, 
— )• 0, and, for every n G N, denote by pn = T„i?„x„. Hence, pk„ x 
and from (3.8) we obtain, for every n G N, 

Pn = TnRnXn <^ Xn - JnBXn £ Pn + InApn 

—{Xn- Pn) - BXn G Apn 
In 
1 

—{Xn-Pn) + Bpn-BXn£{A + B)pn. (3.11) 

In 

Now, since ^ + i? is maximally monotone, from [15, Proposition 20.33], its graph is sequen- 
tially weak-strong closed. Therefore, since Xfc„-pfc„ 0, \\Bpk^-BxkJ\ < x\\xk„-PkJ\ 0, 
7fc„ > £ > 0, pk^ X, we conclude from (3.11) that x G zer(^ + B). Now let us prove that, 
for every i G {1, . . . , m}, fi{x) < 0. Fix i G {1, . . . , m} and, for every n G N, let j„ G N such 
that kn < jn ^ kn + m and i(jn) = i- We deduce from 2;„ — 2;„ — )• and Zn — QnZn — ^ that, 
for every n G N, - Xn\\ = \\QnZn - Xn\\ < \\QnZn - Zn\\ + \\zn - Xn\\ 0. Therefore, 

(Vn G N) \\xj^ — Xk„\\ < — xi\\ < in max \\xi+i — Xi\\ — )• (3-12) 

and hence it follows from Zj^ — Xj^^ — >• and Xk„ x that Zj^ x. Note that, from (3.10) 
and (3.4) we have, for some Uj^ G dfi{zj^), 



(VnGN) Q,„z,„-z,„ = <; ||u,J|2 ' - (3.13) 

0, otherwise, 

and, since ~ -^jnll ~^ 0? deduce that max{0, /j(zj,J} — )■ 0. Thus, it follows from 

Zj^ X that fi{x) < lim fi{zj„) < lim max{0, /^(zj^)} = 0, and hence x G lev<o/j. We 
conclude that x £ Z and the result follows from Theorem 2.2. □ 

Remark 3.6 Let us consider the particular case of Theorem 3.5 obtained when ei^„ = 62, n = 
es^n = 0, m = 1, and /i = dc, where C C ?^ is a nonempty closed convex set. Then, since 
Gi = Pc, Algorithm 3.3 reduces to the method proposed in [7]. Moreover, since S = C, 
note that the assumption ranGi C domi? is equivalent to S C domB, which was already 
assumed in Problem 3.2. 



if/i(^in) > 0; 
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4 Equilibrium problems with convex constraints 



We consider the problem 

find xEC such that (Vy G C) F{x,y)>0, (4.1) 
where C and F satisfy the following assumption. 

Assumption 4.1 C is a nonempty closed convex subset of Ti and F: — >■ M satisfies the 
following. 

(i) (Vx G C) F{x, x) = 0. 

(ii) (V(x, y) G C^) F{x, y) + F{y, x) < 0. 

(iii) For every x in C, F{x, ■): C — )• M is lower semicontinuous and convex. 

(iv) (V(x,y,z) G C^) Wn F{{l-e)x + ez,y) <F{x,y). 

We are interested in solving a more general problem than (4.1), which involves a finite or 
a countable infinite number of convex constraints. It will be presented after the following 
preliminaries. 

Notation 4.2 The resolvent of F: — t- M is the set valued operator 

Jp-.n : x^ {z £ C \ i^y € C) F{z,y) + {z - x \ y - z) > O} (4.2) 

and, for every 5 G ]0, +oo[, the (5-resolvent of F: — M is the set valued operator 

Jl^:n^2^: x^ {zeC \ i^yeC) F{z,y) + {z - x \ y - z) > -6} . (4.3) 



Lemma 4.3 Suppose that F: — )• M satisfies Assumption 4-.1. Then the following hold. 

(i) dom Jp = %. 

(ii) Jp is single- valued and firmly nonexpansive. 

(iii) (Vx G ■H){\/6 G ]0, +oo[) Jpx G ^x. 

(iv) (Vx G ?^)(V(5 G ]0,+oo[) J^px C B{Jpx;VS). 

Proof. (i)&(ii): [20, Lemma 2.12]. (iii): This follows from (ii), (4.2), and (4.3). (iv): Fix 
X ^% and 5 G ]0,+oo[, and let w G JpX. We deduce from (4.2) and (4.3) that F{Jpx,w) + 
{Jpx — X \ w — Jpx) > and F{w, Jpx) + {w — x \ Jpx — w) > —6, respectively. Adding 
both inequalities we obtain F{w,Jpx) + F{Jpx,w) — \\Jpx — > —6. Hence, it follows 
from Assumption (ii) that \\Jpx — < 5, which yields the result. □ 
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Problem 4.4 Let F be a function satisfying Assumption 4.1. Let (S'j)jg/ be a countable 
(finite or countable infinite) family of closed convex subsets of H such that S = Cii^jSi ^ 0. 
Let B : domi? C H ^ Ti he a monotone and x~Lipschitz continuous operator such that 
C C domB, and suppose that 

IJ^i C intdom^. (4.4) 

The problem is to 

find xe S such that (Vy G C) F{x, y) + {Bx \ y-x)>d. (4.5) 

Problem 4.4 models a wide variety of problems including complementarity problems, opti- 
mization problems, feasibility problems, Nash equilibrium problems, variational inequalities, 
and fixed point problems [9, 20, 21, 22, 23, 24]. 

In the literature, there exist some splitting algorithms for solving the equilibrium problem 

find X E C such that (Vy G C) Fi{x,y) + F2{x,y) > 0, (4.6) 

where Fi and F2 satisfy Assumption 4.1. These methods take advantage of the properties of 
Fi and F2 separately. For instance, sequential and parallel splitting algorithms are proposed 
in [25], where the resolvents Jp-^ and Jp^ are used. The ergodic convergence to a solution to 
(4.6) is established without additional assumptions. However, when Fi = F and F2 : {x, y) 1— t- 
{Bx \ y — x) we have = Jb = (Id +B)^^ [20, Lemma 2.15(i)], which is often difficult to 
compute, even in the linear case. Moreover, the ergodic method proposed in [25] involves 
vanishing parameters that leads to numerical instabilities, which make it of limited use in 
applications. In [20, 26] a different approach is developed to overcome this disadvantage 
when B is cocoercive. In their methods, the operator B is computed explicitly and the 
weakly convergence to a solution to (4.5) when 5 = C is demonstrated. 

In this section we propose the following non-ergodic algorithm for solving the general case 
considered in Problem 4.4. This approach can deal with errors in the computations of the 
operators involved. The convergence of the proposed method is a consequence of Theorem 2.2. 



Algorithm 4.5 Let {In)nm be a sequence of finite subsets of /, let (ei,ra)nGN and (e2,n)nGN 
be sequences in Ti such that X^neN < +°'^ Z^nGN ll^2,n|| < +00, and let {6n)neN 

a sequence in ]0, +oo[ such that X^^g^ < +00. Let e G ]0, l/(x + l)[, let (7n)neN 
be a sequence in [e, (1 — e)/x], let U„gN{wi,„}ig/„ C [e, 1] be such that, for every n G N, 



«G/, 



1, let xq G domB, and let {xn)n£N be the sequence generated by the following 



routine. 



(Vn G 



yn = Xn- Jn{BXn + ei,„) 

rn = Qn- ln{Bqn + e2,n) 
— Xn 2/n ~l~ fn 



(4.7) 



Remark 4.6 In Algorithm 4.5, the sequences (ei,„)neN and (e2,n)neN represent errors in the 
computation of the operator B. On the other hand, it follows from (4.7) and (4.3) that, for 
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every n G N, is a solution to 

find q^C such that (Vy G C) F{q, y) + {y - Vn \ V - q) > Sn- (4.8) 

Thus, we obtain from (4.2) that Qn can be interpreted as an approximate computation of the 
resolvent J^^ryn- 

Proposition 4.7 Suppose that there exist strictly positive integers {Mi)i<zi and N such that 

n+Mi-l 

(V(i,n) G / X N) ie |J 4 and 1 < card/„ < TV, (4.9) 

k=n 

and that Problem 4-4 admits at least one solution. Then Algorithm 4-5 generates an infinite 
orbit {xn)n£'N which converges weakly to a solution to Problem 4-4- 

Proof. First, let us prove that Problem 4.4 is a particular case of Problem 1.1. Set 

(Vn G N) /3„ = 7„x, T„ = J^^f, and i?„ = Id -7„S. (4.10) 

Note that {/3n)neN is a sequence in ]0, 1 — e] and, for every n G N, r„ is firmly nonexpan- 
sive [20, Lemma 2.12] and Id —Rn = 7n-B is /3„-Lipschitz-continuous and monotone. Hence, 
it follows from [8, Theorem 1] that the operators {Rn)neN are pseudo contractive. In addi- 
tion, we deduce from (4.2) and (4.10) that (Vn G N) x G Fix TnRn <^ (Vn G N)(Vy G 
C) 'jnPix, y) + {x - RnX \ y - x) >0 <^ (^y £ C) F{x, y) + {Bx \ y - x) > 0. Altogether, 
we deduce that Problem 4.4 is a particular case of Problem 1.1 and 

Z = 5 n Pi Fix TnRn = Sr\{xeC\ i^yeC) F{x, y) + {Bx\y-x)>0} ^ 0. (4.11) 

nGN 

Now let us show that Algorithm 4.5 is deduced from Algorithm 2.1. Set 
(Vn G N) < 

^ Cn = -lne-2-. 

Then, since sup„gpj7„ < we have I]„gNl|a„|| < +oo and I]„gNl|c„|| < +oo. 

Moreover, it follows from (4.7) and Lemma (iv) that X^neN ^ +00, and, for ev- 

ery X ^ TL and n G N, QnX is the projection onto the closed affine half-space Hn{x) = 
{z ^ H I {z — QnX I x — Qnx) < O}, which satisfies S C Rjg/^S'j = Fix Qn C Hn{x) [9, 
Proposition 2.4]. On the other hand, we have xq G domi? and it follows from (4.4) and the 
convexity of intdomi? [27, Theorem 27.1] that 

(Vn G N) ran I ^ uJi,nPs, ] C conv I |J 5, | C conv I |J 5*^ ) C intdomS. (4.13) 

\i£i„ / \ie/„ / Kiel / 

Hence, we conclude from (4.7) that, for every n G N\{0}, Xn G intdomi?. Moreover, for every 
n G N, gri G C C domB. Altogether, from (4.10) and (4.12), we deduce that Algorithm 4.5 
is a particular case of Algorithm 2.1, which generates an infinite orbit {xn)nen- 



an = -7nei,n r ^ _ ^ 

bn = Qn - Jy^pyn and < " „ „ (4.12) 
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Finally, let us show that (2.2) holds. Suppose that Xk„ x, Xn—TnRnXn 0, x„ — 0, 
Zn — QnZn 0, and, for every n G N, denote by pn = TnRnXn- Hence, pk„ x and it 
follows from (4.10) and (4.2) that, for every n G N, 

Pn = TnRnXn ^ (Vz G C) F{pn, z) + — {pn - Xn \ Z - Pn) + {BXn \ Z - Pn) > 

In 

<^ (Vz G C) F(p„,z)H {Pn-Xn\z-Pn) 

In 

+ {BXn - Bpn \ Z-Pn) + {Bpn \ Z - Pn) > 

4^ (Vz G C) G{pn, z) + — {pn - Xn \ Z - Pn) + {BXn - Bpn \ Z - Pn) > 0, 

In 

(4.14) 



where 



G: ^ M: {x,y) ^ F{x,y) + {Bx\y-x) (4.15) 



satisfies the Assumption 4.1 [20, Lemma 2.15(i)]. Moreover, since inf„gN > 0, Xk„ —Pk„ — ^ 
0, and {pkJneN is bounded, we have (Vz G C) {pk„ - Xk„ \ ^ -Pfc„>/7fc„ 0, and from the 
Lipschitz-continuity of B we obtain (Vz G C) {Bxk„ — Bpi^^ \ z — pk„) — ?■ 0. Hence, we 
deduce from pk„ x, Assumption 4.1(iii), Assumption 4.1(ii), and (4.14) that 

(VzGC) G{z,x)<ljmG{z,pkJ 

< lim-G(pfc„,z) 

< lim— (Pfen - Xk„ I z-pkj + {Bxk„ - Bpk„ \ z - pkj 

= 0. (4.16) 

Now let e G ]0, 1] and y £ G. By convexity of G we have Xg = (1 — e)x + ey £ G. Thus, 
Assumption 4.1(i), Assumption 4.1(iii), and (4.16) with z = x^ yield 

= G(x„ xe) < (1 - e)G{x„x) + eG{xe,y) < eG{xe,y), (4.17) 

whence G{xc,y) > 0. In view of Assumption 4.1(iv), we conclude that G{x,y) > 
lim£^o+ G{x£,y) > 0, which yields 

(Vy G C) G{x,y)=F{x,y) + {Bx\y-x)>0. (4.18) 

Now, let us prove that x £ S. Since z„ — — )• and Zn — QnZn — ^ 0, (4.12) yields 

(Vn G N) \\Xn+l - Xn\\ = WQnZn - Xn\\ < WQnZn - Zn\\ + \\Zn - Xn\\ ^ 0. (4.19) 

Now, fix i G /. In view of (4.9), there exists a sequence (jn)neN in such that, for every 
n £ N, kn < jn ^ kn + Mi — 1 and i £ Ij^. For every n G N, it follows from (4.19) that 

\\xj„-XkJ\< y2 \\xe+i - Xi\\ < {Mi - 1) max \\xi+i - xe\\ ^ 0. (4.20) 

<e<k„+AU-2 

Thus, we deduce from xt. ^ x and Zn — x,- — )• that z,- ^ x. On the other hand, let 
z £ S and n G N. Since, for every I £ Ij^, Ps^z = z, and Id —Psi is firmly nonexpansive. 
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from (4.7) and (4.12) we have 

\\Ps,Zj^ -Zjjf < max \\Ps,Zj^ - zjj^ 

< 7 E ^^.i"^^ - I (Id -Ps,)z - (Id -Psjz,^) 



z - Zj„ 



< ^\\z-zjj WQj^Zj^ - z,J\. (4.21) 

Hence, since {zj^)nm is a bounded sequence and Qj„Zj^ — Zj^ — )■ 0, we deduce that PsiZj„ — 
Zj^ — )• 0. The maximally monotonicity of Id —Psi yields that its graph is sequentially weakly- 
strongly closed, and since Zj^ x, we conclude that x = PsiX G Si. Altogether, from (4.18) 
and (4.11) we deduce that x G and the result follows from Theorem 2.2. □ 
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